The dimer-RVB State of the Four-Leg Heisenberg Ladder: 
Interference among Resonances 
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We study the ground state of the 4-leg spin ladder using a dimer-RVB ansatz and the Lanczos 
method. Besides the well known resonance mechanism between valence bond configurations we find 
novel interference effects among nearby resonances. 
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Doped and undoped ladders have focused a lot of at- 
tention in recent years due to the existence of materi- 
als with that structure, some of them are close relatives 
of the high-Tc oxides as the series Sr2ne-2Cu2nt04ni~2 
where is the number of legs forming the ladder [Q. 
The undoped Heisenberg spin ladders with ni even are 
known to be spin liquids with a spin gap and exponen- 
tial decaying correlation functions ||] . The ground state 
(GS) of these low dimensional systems is given by a short- 
range RVB ansatz where the topological spin defects are 
confined j^] . The RVB picture is supported by mean field 
[|, DMRG ||,|, Quantum Montecarlo [§ and Lanczos 
[7| results concerning ladders with = 2,4 legs and 
variational ansatzs (RVA) for the 2-leg ladder |||,^ . The 
purpose of this letter is to apply the RVA method to the 
4-leg spin ladder with the aim of studying in more de- 
tail the structure of the short-range RVB state. In the 
2-leg ladder case the basic mechanism that lowers the GS 
energy is the resonance between two nearest neighbour 
valence bonds . The simplest short-range RVB ansatz 
is given by a dimer-RVB state with a single variational 
parameter u which gives the amplitude of the resonance 
[^,||. In the 4-leg ladder case we shall study a dimer- 
RVB ansatz where the resonance may occur among any 
possible pair of nearest neigbour parallel bonds. The 
new phenomena we shall investigate in this letter is the 
"interference" between couples of resonating bonds. We 
mean by interference the influence that a pair of resonat- 
ing bonds exerts on another pair of near by resonating 
bonds. In the standard RVB ansatz of Liang et al. 
the RVB amplitudes have a factorized form which cannot 
describe this interference effect. 

The Hamiltonian of the 4-leg spin ladder is given by. 



fl" = ^EtiEL"i'Sa(n)-S,(n + l) 



En=i [-/'(SiH • S2(n) + S3(n) • S4(n)) 
+J"S2(n) • S3(n) + J"'Si(n) • 84(71)] 



(1) 



where Sa{n) is the spin 1/2 operator at the a — 1, ... ,4 
leg and n = 1, . . . ,N rung. We shall consider the cases 



of periodic or closed boundary conditions (BC) along the 
rungs, i.e. J' = J" = J'" and open BC's along the rungs, 
i.e. J' = J", J'" = 0. Setting J" = J'" = we recover 
two decoupled 2-leg ladder Hamiltonians. If J = 0, the 
exact GS of (0) is given by the coherent superposition of 
the GS of every rung which can be written as 



|rung) = 12 34 + uq 14 32 

1 J' = J" = J'" 
uo = < 0.366 J' = J", J'" = 



(2) 







J" = J" 







where a6= (| t)a| i)h-| l)a\ T)b)/A/2 denotes the valence 
bond state between the sites on the legs a and b of the 
^th j.^Qg^ jj^ f^gg^ k^qI^ Q^y^^ depict the valence 

bond states (||). 
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FIG. 1. Graphical representation of the local configura- 
tions that make up the dimer-RVB ansatz for the 4-leg ladder. 
Every line connecting two sites a and b corresponds to the va- 
lence bond state ab defined in the text. The site a belongs to 
the even sublattice while the site b belongs to the odd one. 

Switching on the intraleg coupling J any pair of rung- 
bonds will start to resonate with a pair of leg-bonds as 
in figs l(c,d,e,f). There are 4 types of "elementary" res- 
onances involving two consecutive rungs n and n+1 and 
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two legs i and j, which we denote as (12), (34), (14) and 
(23). We associate an amphtude Uij to every of these 
resonances. There is also a state with 4 leg-bonds on two 
consecutive rungs, which we denote as (1234), and give it 
an amplitude M1234 (see fig. 1(g))- Finally, we may have 
two resonances (ij) and (kl) sharing a common rung as in 
figs. l(h,i), which we denote as (12, 34) and (14, 23), and 
give them amplitudes ^12,34 — ^34,12 and 1/14,23 — ^23,14 
respectively. Let us suppose that a pair of resonating 
bonds is not influenced by its environement. This would 
imply the following factorization Uij^ki = UijUu, which 
as we shall see below never happens. Figs. 1 display all 
the local configurations that should be combined in all 
possible manners to produce a dimer-RVB ansatz. This 
seems to be a formidable problem if we try to solve it 
with standard combinatorial methods. However, as in 
the 2-leg ladder case j|] , the dimer-RVB state of the 4- leg 
ladder can be generated by the set of recurrence relations 
(RR) given in fig. 2. 



graph of fig. 4. The vertices of the graph denote the 
configurations A = (0), (12), (34), (14), (23), (1234) while 
a link between the vertices A and A' indicates that these 
two configurations may appear consecutively in an al- 
lowed dimer state. The amplitudes of the dimcr states 
are associated to the links of the graph. 



(12) (0) (34) (12) (0) (23) 

FIG. 3. A dimer state constructed with the RR's given in 
fig. 2. The dotted lines represents the cuts described in the 
text. 
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FIG. 2. Recurrence relations that generate iteratively the 
dimer-RVB state of the 4-leg ladder. The empty box repre- 
sents the singlet state jA'') of the ladder with N rungs. A 
box with two dots on the legs i and j represents the state 
\N,ai, cTj), where and (Jj are free spins which form valence 
bonds with nearest neighbour spins located in the same legs. 
We give explicitely the RR of the state |A'^, o"i, (T2). The RR's 
of the other two-dotted states are similar. The last RR is that 
of the 4-dotted state \N, 0^,02, cry, ,04). 

Fig. 3 shows a state generated by these RR's. One can 
characterize a dimer state with N rungs by the collection 
of legs that one cuts between two consecutive rungs. If no 
legs are cut we write (0), cutting the legs i and j we write 
{ij) and cutting four legs we write (1234). With these 
notations the state of fig.3. reads (12)(0)(34)(12)(0)(23) 
and has an amplitude ui2Ui2,34Ui2W23. 

It is important to realize that not all the states of 
the form A12A23A34 . . . Ajv-i,jv (where A„,„+i denotes 
the set of legs cut between the rungs n and n -f 1) are 
allowed. For example, after the configuration (12) one 
can only have either (0) or (34), or after (1234) only (0) 
may follows. These selection rules are summarized in the 
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FIG. 4. Graph that encodes the dimer configurations of the 
4-leg ladder. The vertices are labelled by the legs cut between 
two consecutive rungs. A link between two vertices represents 
cuts that share a common rung. Every link is associated with 
a variational parameter of the RR's. The link connecting (0) 
to itself means that the middle rung between the two cuts is 
a singlet which may be either 12 34, with amplitude 1 or 
14 23, with amplitude no. 

The RR's of fig. 2 generate all the dimer states of a 4- 
leg ladder with periodic BC's along the rungs, and their 
number grows exponentially with the number of legs p2[ . 
For the open BC's we should restrict ourselves to dimer 
states with no bonds of length greater than one. However 
the strong coupling limit J/ J' << 1 forces to include the 
valence bond 14 as in eq.(||). So the distinction between 
closed and open ladders will only appear in the varia- 
tional parameters. The existence of RR's to generate the 
GS ansatz implies that the norm and expectation value 
of the Hamiltonian (|^) also satisfies RR's, which can be 
iterated to give the energy of the ansatz {H) n for any 
number of rungs A'^. The set of variational parameters 
ux is obtained by minimization of {H)]sf. This method 
is similar to the matrix product ansatz of references |l^] 
but differs in that the states kept are non-orthogonal as 
corresponds to a RVB ansatz. Let us next present our 
results. 
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FIG. 5. GS energy per site in units of J' of the 4-leg 
dimer-RVB state with closed and open BC's and the 2-leg 
ladder in the range < J /J' < 1. We also plot the exact 
GS energies obtained by extrapolating Lanzcos results to the 
thermodynamic limit with ladders of sizes 4 x n(n — 4, 5, 6, 7). 

In fig. 5 we plot the GS energy per site obtained with 
our variational ansatz and the Lanczos method in the 
range of couplings < J/ J' < 1. We also include for 
comparison the GS energy per site of the 2-leg ladder. 
The GS energies are very close to the exact result in the 
strong coupling region < J/ J' < 0.3, but they get worse 
for larger couplings. This is natural since configurations 
with longer bonds are expected to become more impor- 
tant in the weak interleg coupling regime. The closed- 
rung ladder has a much lesser GS energy per site than 
the open one. This is mainly due to the resonance (|^) 
between the two bonds along the rungs. The GS curves 
for open and closed ladders in fig. 5 can be fitted with 
the formula, 

Eo{N)/imj') = -eo - ei(J/J')2 - e2(J/J')' (3) 

, . _ J 0.5, 0.15, -0.005, closed 

(eo,ei,e2j - I 0.23, -0.05 open 

where eo is the energy per site of a single rung. Eq.(|^) 
agrees with perturbation theory up to second order. 

Let us consider now the behaviour of the variational 
parameters. In the closed-rung case the choice of cou- 
plings J' = J" — J'" implies the existence of a rotational 
symmetry among the legs which leaves only 4 indepen- 
dent variational parameters given by, 

UO, Ui2 = Uij, Z;i2,34 = , U1234 (4) 

In fig. 6 we plot these parameters in the domain < 
J / J' < 1. Let us comment on these results. 

• 1*0 takes the constant value 1, which coincides with 
the exact J = result (||) . This implies the absence 
of interference between rung and leg resonance. 
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FIG. 6. Variational parameters M) of the closed BC 
dimer-RVB in the range < J/ J' < 1. We include for com- 
parison the value of the variational parameter u of the 2-leg 
ladder §. 

• U12 is greater than its 2-leg analog u For the 
isotropic case one gets U12 = 1.58 while u = 1.18 

. Simple resonance is enhanced in the 4-leg lad- 
der. 

• ^12,34 is almost constant and less than one indi- 
cating destructive interference between resonances 
shearing a common rung as in fig. 1(c). 

• displays an unexpected behaviour since it first 
becomes negative for small values of J/ J', reaches a 
minimum and starts to grow becoming positive for 
J/ J' > 0.6. This peculiar behaviour of U1234 is a 
sign of destructive interference between resonances 
sharing two rungs. 

In the case of open ladders, J" = J', J'" = 0, one is 
left with 7 independent variational parameters given by, 

Wo, U12 = M34, Ui4, W23, V^j,kl = , ""1234 (5) 

In fig. 7 we plot the values of these parameters in the 
range < J/ J' < 1. Some features that we encounter in 
fig. 7 have already appeared in the periodic case. 

• Uq stays almost constant with a value close to the 
exact J = result (Q). 

• Ui2 and ui4 are quite similar, but U23 is much 
smaller. So bonds do not like to resonate in the 
middle of the ladder. This is due to lost of energy 
induced by the existence of the long bond 14. 

• ^^12, 34 is lower than 1, as in the periodic case, but 
■f^i4,23 is much greater than 1, which is again due 
to the smallness of W23. For graphical purposes we 
plot in fig. 7 the inverse of i;i4,23. 
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FIG. 7. Variational parameters (13) of the open BC 
dimer-RVB in the range < J/ J' < 1. 



• M1234 is also suppressed but in a smaller amount 
than in the periodic ladder. 

We have also computed the spin correlation length ^ 
from the exponential decaying behaviour of the spin-spin 
correlator. For the isotropic case we get ^ = 0.81 for 
the closed ladder and ^ = 0.92 for the open one. These 
results show that the rung configurations are more im- 
portant for the closed ladder than for the open one, which 
is in agreement with the values taken by the variational 
parameters. The DMRG method yields ^ = 5 10 g for 
open ladders, while QMC method yields ^ = 7.1 (closed) 
and ^ = 10.3 (open) As expected the dimer-RVB 
ansatz gives a much shorter correlation length but it re- 
produces the fact that ^closed < 'Copon- 

We have also studied the case when the Hamiltonian 
(|l|) becomes that of two decoupled 2-leg ladders, i.e. 
J" = J'" = 0. Curiously enough our ansatz yields a GS 
with bonds connecting the two ladders. The GS energy 
so obtained is a bit lower than the one of two uncoupled 
2-leg ladders and the correlation length ~ 0.97 is larger 
than in the uncoupled case ^ = 0.79 [||. 

In summary we have shown in this letter that the 
dimer-RVB ansatz gives a correct qualitative picture of 
the short-range RVB state of the 4-leg ladder. We have 
found interesting interference effects between resonating 
valence bond configurations which should probably carry 
over more realistic ansatzs which must include longer va- 
lence bonds. The next step is to generalize our methods 
to the doped 4-leg ladders where one can study the phe- 
nomena of phase-separation and stripe formation P,|l4|. 
Previous applications of the RVA method to the 2-leg 
t — J ladder , the necklace t — J ladder |l^] and the 
2-leg Hubbard model , suggest that this goal is worth 
pursuing. 
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